Vortex qubit based on an annular Josephson junction containing a microshort 
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We report theoretical and experimental work on the development of a vortex qubit based on a 
microshort in an annular Josephson junction. The microshort creates a potential barrier for the 
vortex, which produces a double-well potential under the application of an in-plane magnetic field; 
The field strength tunes the barrier height. A one-dimensional model for this system is presented, 
from which we calculate the vortex depinning current and attempt frequency as well as the interwell 
coupling. Implementation of an effective microshort is achieved via a section of insulating barrier 
that is locally wider in the junction plane. Using a junction with this geometry we demonstrate 
classical state preparation and readout. The vortex is prepared in a given potential well by sending 
a series of "shaker" bias current pulses through the junction. Readout is accomplished by measuring 
the vortex depinning current. 
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I. INTRODUCTION 

Superconducting qubits based on Josephson junctions 
arc one of the most promising qubit architectures in 
terms of scalability and ready integration with semicon- 
ductor electronics. Josephson junctions have been suc- 
cessfully utilized to build various types of qubits such as 
charge^ phased and Aid* 3 - qubits. The operation of these 
systems is based on quantum coherence of the charge 
state, the Josephson phase difference, or the magnetic 
flux state, respectively. Elements characteristic of the 
charge and flux qubits are combined in a hybrid qubit 
called quantronium.— 

By contrast with other types of superconducting 
qubits, a vortex qubit is designed to exploit the coher- 
ent superposition of two spatially separated states for a 
Josephson vortex within a long Josephson junction£&£ 
These states correspond to the minima of a double-well 
potential. The landscape of potential energy experienced 
by the vortex along the length of the junction can be 
constructed as desired by spatially varying one or more 
of the following parameters: the junction barrier thick- 
ness and hence the critical current density^ ' 10 ' 11 ! 12 ? 13 
the magnitude of the in-plane magnetic field jl i 14 i 15 the 
curvature of the junction ccnterline^^^ and the width 
of the junctioni 17 i 18 i 19 

In this article, we report theoretical considerations 
and experimental results concerning a suggested vortex 
qubil£ which consists of a long one-dimensional annu- 
lar junction 2 ^ containing a microshort ^ For this type 
of qubit, the first quantitative fit between an analyti- 
cal model and experimental data in the classical regime 
is presented. The investigated system is schematized in 
Fig. [1] In a microshort qubit, competition between re- 
pulsion at the microshort and pinning by an in-plane 
magnetic field^22 creates a double- well potential for the 
vortex&I 2 -. In having a double-well potential, the mi- 
croshort qubit resembles the flux qubit i 23 > 24 An advan- 
tage of qubit basis states being localized in separate wells 
is that intrawell energy relaxation does not cause a false 



readout result. The quantum state of the microshort 
qubit is manipulated by applying a short pulse of mag- 
netic field, which lowers the height of the microshort- 
induccd potential barrier and thus enables coherent os- 
cillation between the basis states. This field plays the 
same role in controlling the tunneling amplitude as the 
external magnetic flux in the early charge qubit ^ 
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FIG. 1: A vortex qubit consisting of a one-dimensional annu- 
lar junction with a microshort formed by a section of tunnel 
barrier of slightly greater in-plane width. Potential wells de- 
velop to the left \L) and right \R) of the microshort when a 
static magnetic field H, represented by the long arrow, is ap- 
plied. A vortex at rest at either minimum is in a classically 
stable state. The magnetic moment of the vortex is indicated 
by a short arrow. (For clarity, the schematic is not drawn to 
scale.) 

We have enhanced a previously outlined design^ of the 
microshort qubit by implementing a microshort compat- 
ible with standard lithographic fabrication processes.— 
At a "lithographic" microshort, the in-plane width of the 
tunnel-barrier is fractionally larger than it is elsewhere in 
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the junction, resulting in a locally enhanced critical cur- 
rent density per unit length. This planar structure aids 
the monolithic integration of microshort qubits in Rapid 
Single Flux Quantum (RSFQ) circuits, as does the junc- 
tion being large enough for lithographic patterning and 
having a compatible critical current density 3c- With 
RSFQ logic as the interface between vortex qubits and 
room temperature electronics, the circuits could readily 
be scaled up to the large numbers needed for useful com- 
putations. We have previously reported^ that the mag- 
nitude of the bias current required to drive a vortex past a 
lithographic microshort was an order of magnitude larger 
than expected. In this article, we explain that the dis- 
crepancy is due to interaction between the lithographic 
microshort and magnetic field oriented transverse to the 
junction plane. 

A number of characteristics of the proposed microshort 
qubit are tunable during experiment. The applied in- 
plane magnetic field controls the height of the poten- 
tial barrier, the separation of the potential minima, 
and the frequency of coherent oscillation of the vor- 
tex. Also the coupling between microshort qubits is ad- 
justable. Josephson vortex qubits placed in a supercon- 
ducting transmission line couple via virtual electromag- 
netic waves excited and absorbed by vorticespi This in- 
direct interaction depends on the tunneling amplitude of 
each qubit. 

Since quantum tunneling of a single vortex out of a 
metastable potential well has been demonstrated^ 2 - the 
next major step in realizing a vortex qubit is observing 
coherent oscillation of a vortex in a double- well potential. 
Using the experimentally investigated system detailed in 
this article, we have observed for the first time a sin- 
gle vortex escaping from a metastable state by tunnel- 
ing through a microshort-induced barrier^ However we 
limit the experimental data presented here to the classi- 
cal regime, focussing on state preparation and readout. 

The theoretical section of this article, Section [Til be- 
gins with a derivation of the one-dimensional vortex po- 
tential for the proposed qubit. An analytical expression 
is obtained for the depinning current of a vortex over a 
microshort-induced potential barrier as a function of bias 
current and in-plane magnetic field strength. Also, the 
attempt frequency of the vortex in the presence of large 
magnetic field is derived for the cases of zero bias cur- 
rent and bias just below the critical current. From the 
attempt frequency and the height of the potential barrier, 
we find the field dependence of the coupling between de- 
generate potential minima. 

Experimental results obtained in the classical regime 
are presented in Sections IIIII and IIV1 We report mea- 
surements of the vortex depinning current as a function 
of magnetic field strength that indicate the presence of 
bistable vortex states. Our readout scheme, based on 
the unique depinning current of each state, is described 
along with results which confirm that the vortex was pre- 
pared in a chosen initial state by means of shaker bias 
current pulses. The final section, Section IIV1 discusses 



how the vortex depinning current over a lithographic- 
microshort-induced potential barrier can be enhanced by 
flux trapped around both junction electrodes or magnetic 
field applied transverse to the junction plane. 



II. VORTEX IN AN ANNULAR MICROSHORT 
JUNCTION - THEORY 

A. One-dimensional model 




H 



FIG. 2: Plan view of the tunnel-barrier layer of an annu- 
lar junction with a lithographic microshort. The enlargement 
shows the effective microshort in more detail, where the bar- 
rier layer is fractionally wider than it is everywhere else in the 
junction. 

A homogeneous long Josephson junction is well de- 
scribed by the one-dimensional sine-Gordon equation for 
the phase difference ip between the order parameters of 
the superconducting electrodes i 27 i 28 



ip — ip — sin tp = f . 



(1) 



The spatial coordinate x (see Fig. [5]) is in units of 
the characteristic length scale, the Josephson length Aj, 
while time t is normalized to the inverse of the Josephson 
plasma frequency top. The perturbation 



f(x) =7 + 



dh R 
dx 



(2) 



incorporates the bias current 7 and the current induced 
by the radial component hp of the normalized external 
magnetic field h = kH/Hq. The bias current 7 is nor- 
malized to the product of the critical current density Jc 
and the junction area. The characteristic field is given by 
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Hq = $o/27r^o<iA,/ where $o is the magnetic flux quan- 
tum, p, is the vacuum permeability, d is the magnetic 
thickness of the junction, and k is the geometric cou- 
pling factor between the magnetic field applied and that 
in the junction. For a linear junction dhu/dx = 0, but 
in the annular case dhu/dx = — (h/r) sin(x/r) where r 
is the junction radius. Dissipative terms have not been 
included in Eq. |TJ because, for the purposes of the ex- 
periments described in Sections IIIII and HVl the damping 
in our system at low temperatures is small enough to 
neglect. 

In the absence of external field and applied bias cur- 
rent, one of the solutions of Eq. {1} is that of the station- 
ary Josephson vortex: 



<P 



4 arctan[exp(x — xq)] 



(3) 



The Josephson vortex is a topological soliton, behaving 
like a particle of normalized rest mass energy too = 8 
and center of mass position xo- Here the unit of energy 
£o = Jc w o^j'&o/'Z'R is the Josephson coupling energy of a 
small junction of area wq Xj, wq being the junction width. 
To find the potential seen by a vortex in an annular 
junction, we begin with the Hamiltonian corresponding 
to Eq. ©: 



Hq 



(\ 9 1 2 \ 

\2 V " + 2^ + 1 ~ ° OS ^ + V ' 



(4) 



Approximating the spatial phase profile <p(x) in Eq. (j4|) 
with Eq. (J3]), that of a vortex, and carrying out the in- 
tegration over the entire length I = 2irr of the junction, 
one obtains the washboard potential:— 



U(xq) = — 2tt[jxq + hsech — cos — 

2r r 



(5) 



The first term U 1 ~ —2ttjxq reflects the driving ac- 
tion of the bias current on the vortex, while the sec- 
ond, U h = — 27r/isech(7r/2r) cos(xo/r), originates in the 
convolution of the vortex spatial magnetic field profile 
dip v /dx = 2sechx with the applied field component Hr: 



jjh 



d(p v 
dx 



(6) 



Next we extend Eq. ([5]) by adding the contribution due 
to the lithographic microshort depicted in Fig. [21 

The characteristic energy scale £q and hence the 
Josephson vortex rest mass energy m are proportional 
to the width of the junction. Thus one expects that a 
short length Z M of broader junction will act as a potential 
barrier to the vortex, just as a microshort formed by de- 
creasing the tunnel barrier thickness would. In the case 
where the lithographic microshort length is smaller than 
the characteristic vortex size, l^ < Xj, it is important to 
consider the spatial distribution of the vortex mass. 

The spatial rest mass profile is found from the sine- 
Gordon Hamiltonian for a static vortex: Eq. (0| with 
/ = 0. The time- independent phase distribution, Eq. ([3]), 



gives rise to equal magnetic field and Josephson coupling 
energy densities, f' 2 /2 and 1 — cos(p respectively. Hence 
the vortex rest mass energy can be written as 



toq 



jldx 



(7) 



where the spatial distribution of the vortex rest mass is 

jl — 4sech 2 (x — x ) . (8) 

To calculate the change of vortex rest mass energy 
caused by altering the junction width, we start with an 
appropriate Hamiltonian 



H 



w{x) /l 



W 



2*' 



1 — cos ip ) dx 



(9) 



neglecting the perturbations due to bias current and ex- 
ternal magnetic field. Note that the energy normaliza- 
tion Eq remains constant since the factor w(x)/wq de- 
scribing the spatial variation of the junction width ap- 
pears explicitly in Eq. ([9]). We only consider small width 
changes, so the phase distribution ip(x, y) = <p(x) remains 
radially independent r^. One now sees that the vortex rest 
mass energy mo in a long junction of variable width con- 
sists of the convolution of the spatial rest mass profile 
belonging to a uniform junction with the local potential 
energy per unit length of the variable width junction: 



TO = /i : 



Wq 



(10) 



The kernel jx ~ sech (x — xo) of this convolution in- 
tegral has a different shape to the magnetic field pro- 
file dip v /dx ~ sech(x — xo) in Eq. ([6]). The rest mass 
and magnetic field profiles are plotted together in Fig. [31 
For the same amplitude, the rest mass profile is nar- 
rower. This means that the potential barrier caused by 
a microshort tends to be more spatially confined than 
a barrier produced by an external field. Consequently, 
quantum tunneling of a Josephson vortex through a 
microshort-induced barrier is expected to be enhanced.— 

Turning to the specific long junction geometry depicted 
in Fig. [21 the width variation is 



w(x) 

Wq 



Wq 



-n(x/y 



(ii) 



where n(x) is the unit rectangle function and Aw = w^~ 
wq is the amount by which the junction width at the 
lithographic microshort, w^, is larger than the junction 
width elsewhere. For short microshorts, l^ ^C Aj, an 
equivalent representation using the Dirac delta function 
is 



w(x) 

Wq 



1 



Aw 

Wq 



lfj,S{x) 



(12) 



Then the potential energy corresponding to the change 
in vortex rest mass energy evaluates to 



{/^(xq) = /j, sech 2 Xq 



(13) 
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(a) 100 



FIG. 3: Comparison of the kernels dip v /dx and jl of the con- 
volution integrals which respectively give rise to the magnetic 
field U h (xo) and microshort U ,J '(xo) contributions to the vor- 
tex potential. The rest mass profile jl ~ sech 2 (a; — xo) (solid 
line) of the vortex is spatially more tightly confined than 
its magnetic field profile dip v /dx ~ sech (a: — xo) (dashed 
line); Therefore the resultant microshort barrier tends to be 
narrower, and the associated vortex quantum tunneling rate 
should be enhanced.— 



where /i = AI^Aw/wq measures the strength of the mi- 



cro-short. Combining Eq. (fl~3|) and Eq. ([5|) gives the vor- 
tex potential for an annular junction containing a mi- 
croshort: 



U(xq) = /isech .To — 2tt( 



7^o 



h sech — cos 
2r 



*)•<"> 



B. Double well potential 



The vortex potential in the absence of bias current, 
Eq. (|14[) with 7 = 0, has two wells for magnetic fields in 
the range < h < 1 where 



hir sech(-7r/2r) 
fj,r 2 



(15) 



Figure 2£i contains plots of this potential for various val- 
ues of the magnetic field. The parameters [i, Aj, and r 
were chosen to reflect the experimentally investigated 
system, which is detailed in Section IIII Al The height 
of the potential barrier shrinks with increasing magnetic 
field strength. In this way, the barrier height can be 
quickly controlled during experiment by, for example, 
passing current through an appropriately oriented mi- 
crostrip underneath the junction. The field at which 
the barrier disappears completely is found by calculat- 
ing when the curvature U"(xq) of the potential becomes 
positive at xq = 0. The left \L) and right \R) wells of the 
potential constitute stable classical states for the vortex. 
The states are distinguished by their critical, or de- 
pinning, current. This is the amount of bias current 




FIG. 4: The vortex potential U(x ) given by Eq. (fT4|) for 
a long annular junction of radius r = 3.7 containing a mi- 
croshort of strength /i = 0.09 . (a) The potential is plotted 
for four magnetic fields in the range ^ h ^ 1 at zero bias cur- 
rent. The double- well potential is tunable: As the magnetic 
field strength increases, both the height of the microshort- 
generated potential barrier and the physical separation be- 
tween the stable vortex states \L) and \R) decrease. A black 
circle marks the location xc2 at which an inflection point 
forms under the critical bias yc2- (b) The potential at h — 0.8 
is shown for four bias currents in the range 0^7^ 702- A 
white circle indicates the location of the left minimum, which 
becomes an inflection point at 7 = 702 • 



needed for the vortex to overcome the pinning poten- 
tial. We denote the depinning currents for the field- and 
microshort-induccd barriers as jci and 7C2, respectively. 
As shown in Fig.[3J) for field h = 0.8, a larger bias current 
tilts the potential further. The left potential minimum, 
indicated by a white circle, turns into a horizontal in- 
flection point at the critical bias "fc2- The location at 
which the inflection point develops, Xc2, is represented 
in Fig. 2£i by black circles. The depinned vortex gains 
kinetic energy as it moves through the microshort and 
beyond in the direction of increasing x, and the junction 
voltage becomes nonzero. The vortex is metastable at 
bias currents slightly below that required to transform 



C Attempt frequency 



its potential well into an inflection point. This is due to 
thermal fluctuations and quantum tunneling, which re- 
duce the observable critical current. At a large enough 
magnetic field h, a vortex which has depinned from the 
left \L) well will be retrapped by the field-induced po- 
tential barrier before the switch to the resistive state can 
be observed. This can be avoided during readout of the 
qubit state by decreasing the field before the depinning 
current is measured. 

We now consider the critical current lci(h) associated 
with a vortex escaping over the potential barrier induced 
by the magnetic field. Increasing positive bias current 
moves the right minimum \R) and the vortex farther in 
the positive x direction. At the critical bias, the vor- 
tex reaches the location xci = 7rr/2, where its magnetic 
moment is perpendicular to the applied field and the po- 
tential has a saddle point. Here the vortex depins and the 
bias accelerates it. In the case of a long annular junction 
with a weak microshort strength /j, <C 1, it is reasonable 
to neglect the influence of the microshort on this critical 
current. Therefore^ 



lci(h) 



- — - sech — 
r 2r 



(16) 



Next we derive the critical current 102(h) due to the 
microshort. In this case, we assume xq ~ and replace 
the cosine term in Eq. ([T4"[) with its Maclaurin series up 
to O(xq), neglecting an additive constant: 



U(xq) — n (sech xq 



■fuel) 



2ttjxq . 



(17) 



At the critical bias 7c2, the vortex depins at the loca- 
tion xci where the potential has a saddle point: 



U'(x C2 ) = and 
U"(x C2 )=0. 



(18) 
(19) 



The condition given by Eq. (TH)|) leads to a quadratic 
equation in sech xc2 with the solution 



sech xci 



xci 



1+ Vl + 3h 



3 

arcsech 



n(h) 



(20) 
(21) 



The minus sign in Eq. (|21|) reflects the case of a vortex 
which moves in the positive x direction as it escapes over 
the microshort-induced potential barrier. A larger ap- 
plied field in the range < h < 1 results in a depinning 
location xci which is nearer to the microshort center 
at x = 0. After substituting Eq. (JUl) into Eq. fH]), one 
obtains the field-dependent critical current for a vortex 
trapped by a microshort: 



/' 



lci{h) = — (tiVI - n - h arcsech y/rj) 



(22) 



C. Attempt frequency 

A vortex in a potential well oscillates about its aver- 
age position with frequency u>q, known as the small am- 
plitude oscillation frequency. This is the frequency with 



which the vortex attempts to escape from the well; In a 
junction with moderate damping it is proportional to the 
rate of thermal escape. [See Eq. ([42]) .] In the quantum 
regime, the energy hu>o separates the first excited state 
of the vortex from its ground state within a single well. 
For a symmetric double-well potential, the frequency of 
coherent oscillation of the vortex between the \L) and \R) 
states, Ao, depends exponentially on the small oscillation 
frequency oj [as seen from Eqs. ([36]) and ([37]) ]. In this 
section, we derive analytical expressions for the attempt 
frequency coq at large fields h < 1. 

The frequency of attempted vortex escape varies as the 
square root of the curvature at the potential minimum: 



w 



f £/"(zi) 

m 



(23) 



The location of the potential minimum is found by solv- 
ing 



U'{xi) = 



(24) 



using the approximation x ~ tanh x + tanh x/3, which is 
valid for small X{ and therefore large normalized fields h. 
The approximation leads to an analytically solvable cubic 
equation in tanna^: 



tanh Xi + t-2 tanh x% = tj, 
where, for clarity, we define 



Tl(h) 



1 



l + h/3 



and T 3 (h) 



7T) 



fi(l + h/3) 



(25) 



(26) 



The discriminant of the cubic equation, D = (T2/3) 3 + 
(r 3 /2) 2 , equals zero at the critical bias 702, where all 
roots are real and two are equal. From this, one obtains 
an approximation to Eq. 



102(h) 



2 M (1 - ~hf' 2 
3V3tt(1 + h) 1 / 2 



(27) 



whose relative error is under 0.7% over the range of mag- 
netic fields for which the double-well potential exists. 
Defining the normalized bias current as 7 = l /102(h) 
with ~iC2(h) given by Eq. ([27[) enables the solutions 
to Eq. ([25]) to be written as 



, -Ti /arccos7 2n7r\ 

/ -± cos f — '- + — J . (28) 

The constant n takes on values from the set {0, 1, 2} for 
a double-well potential where the microshort is centered 
at x = 0. The set elements correspond respectively to the 
right potential minimum \R) at X3, the left minimum \L) 
at xi, and the potential maximum at X2- 

An analytical expression for the attempt frequency 
is found in the case of zero bias current; Substitut- 
ing Eq. (gBD with 7 = and n = 1 into Eq. Q leads 
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to 



U! 



yjlih (1 - ft) (33 - h) 



2 (3 + h) 



(29) 



For large fields h < 1, near the field at which the 
microshort-induced potential barrier disappears, a series 
expansion of the square of Eq. (f29|) about h = 1 reveals 
that the attempt frequency behaves as 



w 



v7"W 2 



(30) 



where e/i = 1 — /i. 

The effect of bias currents just below the critical cur- 
rent, 7 < f , on the attempt frequency is of interest when 
interpreting microwave spectroscopy data for a vortex in 
a metastable state. Such data provide information on 
the shape of the potential well. For large magnetic fields, 
h < 1, and large bias currents the curvature of the po- 
tential minimum can be described by 



£/"(zi) ~ 4 M e /l y / 2^/3 



(31) 



where &y = 1 



7. On inserting Eq. J5IJ into Eq. (|2"3"j) . 
one finds that the attempt frequency for a vortex pinned 
by a microshort in a metastable potential well depends 
on the bias current as 



where U(xq) is given by Eq. (fT7|) . 

At zero bias current, the tunnel splitting ftA depends 
on the instanton action S as2£ 



^H exp H 



where 



S 



16AC/ 

3cjn 



(36) 



(37) 



and AC/ denotes the height of the potential barrier. 
From Eq. (17|) and Eq. ([28]) with ra = 1, 



AU = ^(tanh 2 x x - hx\) 



3/ze£/8 for h < 1 



(38) 
(39) 



Hence the action varies with large magnetic fields ft as 
S = 2v/2^ 2 , (40) 

and the quantum tunneling rate Ao, in units of up, is 
given by 



/-^f/^cxp 



2V2>f N 



(41) 



(32) 

Just as for a vortex pinned by a microresistor in the 
absence of magnetic field^^ the small oscillation fre- 
quency wo varies with bias 7 as the fourth root of the 
term 1 — 7. 



D. Quantum properties 

To investigate the quantum properties of an annular 
Josephson junction with a lithographic microshort, we 
start with the normalized Euclidean action S, in units 
of£o/ujp 7 at zero temperature: 

/>oc 1*1/2 

S\ip{x, t)] = / dr dx 

Jo J-1/2 

,2 



tf 



1 



S(x) 



2 



+ 1-008^ +ff> ■ (33) 



Here the time variable is transformed as t = it, and the 
vortex coordinate xq(t) is time dependent: 



tp v (x, t) = 4arctan{exp[x — xq(t)]} . 



(34) 



After integration of Eq. (|33|) over the spatial coordi- 
nate x, the effective action has the form 



S[(f(x,T)] 



dr 



m io 



U(x ) 



(35) 



E. Operation as a two-level system 

Coherent oscillation of the vortex between the basis 
states \L) and \R) is controlled by short pulses of the in- 
plane field H . At zero bias current, these states are de- 
generate. Oscillation takes place during a positive pulse, 
when the potential barrier AC/ is lower and the interwell 
coupling Ao is thereby increased. At this field, the po- 
tential minima are close enough to each other that the 
crossover temperature To, below which exponential relax- 
ation gives way to undcrdamped coherent oscillation, is 
greater than the temperature T of the qubiti^S After the 
end of the pulse, the crossover temperature To is lower 
than the qubit temperature T, and the oscillation has 
ceased. Readout should take place before overdamped 
relaxation causes the vortex to transition to the other 
well. At this point, there is enough time to further lower 
the field if it is necessary to reduce the barrier height 
before the depinning current of the \L) state can be mea- 
sured. 

We now consider the feasibility of experiments in the 
quantum regime using realistic trilaycr fabrication pa- 
rameters. We assume a critical current density of Jc = 
1 kA/cm 2 , a junction width of wq = 0.6 /im, and an ad- 
ditional junction area of (w^ — Wo)l)j, = 0.2 (im 2 . This 
corresponds to a Josephson length of Xj « 12 /xm, a mi- 
croshort strength of fj, = 0.1, and a reduced Planck's 
constant of ft, — 3.5 x 10 -3 £o/wp. I n normalized units, 
the magnitude of the reduced Planck's constant ft is in- 
versely proportional to the junction width wq. Hence 
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FIG. 5: Quantum tunneling rate Ao (solid line) and up- 
per limit for the thermal activation rate I\h (dashed line) 
for a junction of critical current density Jc = lkA/cm 2 , 
width wo = 0.6 /im, and microshort strength /j, — 0.1. 



quantum effects are more readily measured in narrower 
junctions^ The length of the junction is taken to be 
30 Aj. 

We estimate an upper limit for the thermal escape 
rate T t h by assuming moderate damping: 



r tll = -exp(-^j 



(42) 



The actual thermal activation rate is expected to be sub- 
stantially slower because the junction will be in the low 
damping regime. Figure [5] aids comparison of the up- 
per limit for the thermal escape rate T t h at temperature 
T — 50 mK with the quantum tunneling rate Ao given 
by Eq. (|4"Tj) . Clearly, thermal activation should be in- 
significant during operation of the junction as a two-level 
system at T = 50 mK. 

The friction coefficient r\ plays a crucial role in the 
dissipative dynamics of a two-level system. It is given by 



1 



77 = 



ujpRC 



(43) 



where C denotes the junction capacitance and R = 50 0, 
the real part of the impedance of the microstrip biasing 
lines feeding the junction. For our chosen junction pa- 
rameters, the classical friction coefficient is 7/ = 2 x 10~ 3 . 
To determine the experimental conditions under which 
coherent oscillation of the vortex between the left \L) 
and right \R) potential minima is observable, we use the 
dimcnsionless system-environment coupling strength 



?7(Aa;) 2 
2irh 



(44) 



where Aa; is the distance between the minima according 
to Eq. l|2"8"|) . The crossover temperature Xq depends on 



FIG. 6: Barrier height AU (solid line); spacing huio (dashed 
line) between the first excited state and the ground state in 
each well; and crossover temperature To (dotted line), be- 
low which exponential relaxation becomes underdamped co- 
herent oscillation. The junction parameters are the same as 
for Fig. El 



the coupling strength a a 



T 



,30 



_fi_Ap_ 



(45) 



Figure [5] displays the crossover temperature Tq, the bar- 
rier height AC/, and the spacing Hujq between the first 
excited state and the ground state in each well. At a 
temperature of T = 50 mK, an applied field of e^ ~ 0.06 
for the quantum operation and e^ ~ 0.14 for the read- 
out should yield a measurable quantum tunneling rate of 
Ao ~ 700 MHz and allow a readout time of up to tens 
of milliseconds. Rectangular field pulses of the required 
height Ah = 0.08 and frequency could be produced using 
a microstrip or the near field of an rf antenna. During 
quantum operation at eh ~ 0.06, the level spacing fitoo 
is one order of magnitude greater than the tunnel split- 
ting ftAo- Thus pulse frequencies comparable to the tun- 
neling rate Ao can be used without excitation of the vor- 
tex to a higher level provided that the rise time of the 
field pulse is longer than 2-k /ojq. 



III. EXPERIMENT 
A. Samples and setup 

The results presented here pertain to two nominally 
identical long Josephson junctions of the geometry de- 
picted in Fig. [21 The junctions were fabricated on differ- 
ent chips within the same run of a standard Nb-AlOz- 
Nb trilayer process^ and had a critical current density 



of Jc = 1.2kA/cm at temperature T = 1.4 K, which 
corresponds to a Josephson length of Xj = 13 /im (after 
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including the effect of the idle region)^ The basic shape 
of the junctions was annular, with a radius of r = 50/im 
and a width of wq = 2.4 /im. A short length l^ = 1.4 /jm 
of each junction was wider, w M = 2.9 (im, forming a 
lithographic microshort of expected strength /x = 0.09. 
The junctions had Lyngby style^ electrode leads, whose 
width slightly exceeds the junction diameter 2r. 

Each sample was situated within a solenoid which pro- 
duced the in-plane magnetic field H . For measurements 
of the vortex depinning current Ic{H) at T = IK, the 
solenoid was attached to the cold finger of a dilution re- 
frigerator. The junction was biased through low pass 
RC filters at the 1 K staged Measurements at T = 1.4 K 
were carried out in a 4 He cryostat. For both setups, room 
temperature feedthrough capacitor filters were used. The 
Josephson vortex was trapped upon cooling through the 
critical temperature of the niobium electrodes. Its depin- 
ning current was evaluated as the product of the known 
bias current ramp rate with the time elapsed between the 
ramp passing through zero current and the appearance 
of voltage across the junction electrodes^ 



B. Bistable states 

Measured data on the dependence of the vortex de- 
pinning current Ic as a function of the external mag- 
netic field strength H are plotted in Fig. [7] together with 
Eqs. flTEJ) and d22J). The I C (H) data consists of lines 
which are labelled according to the source of the poten- 
tial barrier over which the vortex escaped: i) the mag- 
netic field, ii) the lithographic microshort, and Hi) the 
injectors. 35 ' 36 The branches which appear at high field H 
originate from the nuclcation of iv) one or v) two vortex- 
antivortex pairs.— The injectors are leads for local cur- 
rent injection, which provide a means of inserting a vor- 
tex into the junction on demand. They were not used 
in the experimental work described in this article and 
they do not affect the shape of the relevant part of the 
double-well potential, even though the injector branch is 
present at weak positive fields H. This is because the 
current injection leads are located at x = 7rr, as far from 
the microshort as possible. 

The dependence of the depinning current on the 
magnetic-field barrier is typical for long angular junc- 
tions. At zero field the barrier is absent and the junc- 
tion critical current, at which the vortex depins, is there- 
fore minimal. Increasing the magnetic field strength 
leads to a larger vortex depinning current as expected 
from Eq. (|16p . which is plotted as a solid line in Fig. [7] 
In order to graph Eq. (flU)) , the geometrical parameter k 
was calculated from the gradient of the branches labeled i 
while taking into account the self-field of the bias current. 
The maximum critical current was measured indepen- 
dently, and the values of the other junction parameters 
came from the specifications of the fabrication process. 

Linear extrapolations of branch i for positive and nega- 
tive field H do not intersect at zero critical current. This 




FIG. 7: Vortex depinning current Ic for junction #1 as a 
function of applied magnetic field strength H at tempera- 
ture T = 1.4 K. The standard deviation, which is on the order 
of 1 jj,A, is not shown. Sections of the data (black circles) form 
approximately straight lines which originate in the escape of 
the vortex over a potential barrier induced by either i) the 
magnetic field, ii) the microshort, or Hi) the injectors . 35 ' 36 
The depinning currents expected from Eqs. (|16|l and (|22|l for 
the magnetic field and microshort barriers are also plotted 
(solid and dot-dashed lines, respectively). In calculating the 
depinning current from the magnetic field barrier, the self- 
field of the bias current was taken into account. The inset 
displays the supercurrent maxima of the Ic(H) pattern as 
well as branches due to the nucleation of iv) one or v) two 
vortex-antivortex pairs. 



could be due to the self-field of a current induced by a 
small magnetic field Hz oriented perpendicular to the 
junction plane. The induced current circulates along the 
electrode rings. At xq = +nr/2 and xq = —irr/2, where 
the vortex depins from the barrier induced by positive 
and negative field H, respectively, the self-field of the 
circulating current is of opposite polarity. The contribu- 
tion of this self-field to the total in-plane field H shifts 
the measured branches i, for both positive and negative 
field H, along the H axis towards the origin. 

For negative field polarity, the magnetic-field barrier is 
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centered at the microshort, at xq — 0, which results in a 
single-well potential. For positive polarity, the magnetic- 
field barrier, at xq = irr, is located diametrically op- 
posite the microshort barrier; Together they form the 
double-well potential. As seen in Fig. [Jj at temperature 
T = 1.4 K the stable vortex states of the double-well 
potential are distinguishable by means of the depinning 
current over the field interval < H < 5 A/m. Depin- 
ning currents from both of the two branches i and ii are 
also measured for fields near H = —4 A/m although the 
potential at zero bias consists of a single well. In this 
case the junction radius was large enough that, at finite 
bias, a second local potential minimum formed between 
the positions where the respective gradients of the field 
and microshort contributions to the potential, dU h jdx 
and dU^ /dx, were maximal. 

In Fig. [7J the measured depinning current from the mi- 
croshort barrier is several times larger than predicted by 
Eq. ((22)) . which is represented by a dot-dashed line. We 
discuss this issue in Section HVBl where we show how an 
external magnetic field Hz perpendicular to the junction 
plane can result in an effective microshort strength /2 
which grows with increasing field strength Hz- When 
ft 76 //, the dependence of the vortex small oscillation fre- 
quency ojo on the in-plane field h differs from Eqs. (|2"5)) 
and (|30p . This prevents direct comparison of the data 
in Fig. [S] and Fig. [JJ One means of avoiding interaction 
between the microshort and transverse field Hz is sug- 
gested in Section HV CI 



C. State preparation 

To base a quantum bit on this junction design, it is 
necessary to be able to produce a defined initial state. 
We prepared the vortex as desired in the left \L) or 
right \R) potential well by means of a shaker— bias cur- 
rent sequence,— which is depicted in Fig. [8^,. The ampli- 
tude \Ip\ of the preparation pulses is such that a pulse 
frees the vortex only when it is pinned by the magnetic 
field and not when it is pinned by the microshort. The 
vortex remains pinned in the left \L) and right \R) wells 
for pulses of positive and negative polarity, respectively. 
With each successive pulse, the likelihood of locating the 
vortex in the chosen well increases. The number of pulses 
needed to reliably attain the desired state depends on the 
probability of the vortex being retrapped in the chosen 
well. Note that the vortex can also be prepared deter- 
ministically in the left \L) or right \R) well by applying a 
bias pulse of the appropriate polarity at zero field and in- 
creasing the external field afterwards, but this procedure 
is slower due to solenoid inductance. 

After the preparation stage, the bias current is in- 
creased at constant positive rate in order to measure 
the vortex depinning current and identify which poten- 
tial well the resting vortex was located in. The readout 
results^ for thirty preparation pulses of amplitude | Ip \ = 
1.44 mA in junction #2 at temperature T = 1 K are plot- 
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FIG. 8: State preparation by means of a series of shakes^ 
bias current pulses applied to junction #2 at temperature T — 
1 K.— (a) The bias current during state preparation and read- 
out, (b) Positive (white triangles) and negative (grey cir- 
cles) preparation pulses place the vortex in the left \L) and 
right \R) potential wells, respectively. The magnitude of the 
preparation pulses (dashed line) is indicated. 



ted in Fig. [8)3. The white triangles denote depinning cur- 
rents measured after preparing the state with positive 
polarity pulses, and the grey circles, negative polarity 
pulses. 

We observed over tens of trials that thirty preparation 
pulses are enough to prepare a given state at nearly all 
external magnetic field strengths. An exception occurred 
near 10 A/m where the vortex could not be prepared in 
the left well \L) due to the low probability for it to be 
retrapped there. We found that the probability for the 
vortex to be retrapped in a particular well is tempera- 
ture and field dependent. The temperature dependence 
of the retrapping probability may be caused by the damp- 
ing changing with temperature. Initialization at a field H 
where the probability of retrapping the vortex in the de- 
sired well is high offers the advantage that fewer prepa- 
ration pulses are required. We expect that the shaker 
initialization process can be used at most fields where 
the depinning currents from the \L) or \R) wells are dis- 
tinct as long as the amplitude \Ip\ of the preparation 
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pulses lies between these currents and is greater than 
the depinning current associated with any other poten- 
tial barriers in the junction. Therefore junction #1 could 
be initialized with this procedure over a field range from 
H = A/m to H ~ 5 A/m, a shorter interval than for 
junction #2. 

In Fig. [SJd, the measured depinning current from the 
left well \L), where the vortex is pinned by the mi- 
croshort, is an order of magnitude greater than predicted 
by Eq. (f22|) . Such a large enhancement occurs when 
magnetic flux is trapped around both electrodes. The 
flux-induced supercurrent interacts with the lithographic 
microshort as described in the next section. 



IV. DISCUSSION 



Enhanced vortex pinning due to flux trapped 
around both electrodes 



We have found that the vortex depinning current 
strongly depends on the strength of a magnetic field Hz 
applied perpendicular to the junction plane during cool- 
ing below the superconducting transition temperature 
Tc ~ 9.2 K. After field cooling the transverse field Hz 
was switched off, and then the critical current Ic was 
measured as a function of the in-plane field H. We re- 
peated this process for a range of transverse fields Hz 
and obtained several distinct Ic{H) patterns. Each 
Ic(H) pattern was observed for a particular subinterval 
of transverse field. Reversing the polarity of the applied 
transverse field reflected the Ic(H) patterns about the 
Ic and H axes simultaneously. Further results pertain- 
ing to field cooling annular junctions under perpendicular 
fields are presented elsewhere i 39 ' 40 

Figure [5] shows Ic (H) patterns which were obtained 
after field cooling under transverse fields Hz differing by 
around 0.2 A/m and equivalent to fluxes of $z = —2.8 $0 
(grey) and <&z = — 3.7$o (black) through each electrode 
loop. Whereas the Ic(H) pattern obtained after field 
cooling under $2 = —2.8 $0 is typical of a single trapped 
vortex, the pattern for &z = — 3-7 $0 resembles Fig. [SJd 
in that the vortex depinning current over the microshort- 
induced barrier is enhanced by an order of magnitude. 
The difference between these two Ic(H) patterns origi- 
nates in the amount of flux trapped within both electrode 
rings. Note that the microshort depinning current mea- 
sured after field cooling under <&z = —2.8 $0 does not 
agree with Eq. (f2"2"| . possibly due to the influence of a 
weak background transverse field. 

Field transverse to the junction plane induces a super- 
current /$ which circulates around the electrode loop. 
The fraction of supercurrent in each electrode is in the 
inverse ratio of their inductance per unit length. Due 
to the fabrication process, this ratio changes at the mi- 
croshort boundaries, where not only the active region but 
also the neighboring part of the upper electrode widens as 
illustrated in Fig. [TJ Consequently, some supercurrent e 
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FIG. 9: Vortex depinning current Ic versus in-plane magnetic 
field H at T — 1.4 K, obtained after field cooling junction #1 
under transverse magnetic fields equivalent to —2.8 (grey) and 
—3.7 (black) flux quanta through each electrode loop. The 
perpendicular field was switched off before the Ic(H) pattern 
was recorded. 



passes upwards through the junction barrier on one side 
of the microshort and downwards on the other. The cor- 
responding spatial dependence of the Josephson phase 
ip(x) is equivalent to that of a current dipol o 35 ' 36 or frac- 
tional vortex ; 41 i 42 objects which can increase the vortex 
depinning current over the microshort-induccd potential 
barrier. For a pointlike (L <C A,/) dipole of strength 



dfi, Eq. @ become; 



,36 



dh. 



/(x) = 7+^f + M'0*0- 



(46) 



The contribution of the lithographic microshort to the 
vortex potential is then 



U fl (xo) = /xsech xq — 2« M sechxo 



(47) 



= (ji - 2k m ) sech 2 x - ft^o + °( x o) ■ ( 48 ) 

A negative dipole strength k m , which corresponds to the 
induced supercurrent e circulating in the direction of de- 
creasing x, thus increases the height of the potential bar- 
rier at the lithographic microshort. 

The order-of-magnitude enhancement of the mi- 
croshort depinning current evident in Fig. [SJd is due to 
the dipole current e induced by magnetic field perpendic- 
ular to the junction plane. The field is from flux quanta 
trapped within the electrode rings, possibly in combina- 
tion with a residual background field Hzo- The back- 
ground field for the data presented in Fig. [7] and Fig. [8] is 
likely to be different because they were recorded in dif- 
ferent cryostats. An increase in the depinning current 
due to trapped flux takes on discrete values, whereas the 
increase caused by a background field does not and can 
therefore be much smaller. Another difference in the ef- 
fect of trapped flux and a background field is that they 
induce dipole currents e of opposite polarity. 



A Enhanced vortex pinning due to flux trapped around both electrodes 



11 



B. Enhanced vortex pinning due to magnetic field 
applied transverse to the junction plane 
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FIG. 10: Vortex depinning current Ic as a function of ap- 
plied in-plane magnetic field H and perpendicular field Hz- 
Measurements were recorded at a temperature of T = 1.4 K 
using junction #1. The microshort depinning current at 
H = OA/m depends linearly on the applied perpendicular 
field Hz- Also plotted is a linear fit (solid line) to the mag- 
netic field depinning current for Hz = OA/m. The fitted 
depinning current from the microshort barrier is indicated 
for two cases of misalignment of the solenoid field H : 0.057° 
out of (dotted line) and 35° within (dashed line) the junction 
plane. The inset displays the Ic(H) pattern as a function 
of Hz over a field H interval which includes the supercurrent 
maxima. 

We have investigated the effect on the microshort de- 
pinning current of a uniform external field Hz applied 
transverse to the plane of a junction in the superconduct- 
ing state. The transverse field Hz was generated using a 
current coil located underneath the chip. Note that, since 
the field lines produced by the coil are screened by the 
superconducting electrodes and cannot pass through the 
hole of the electrode loop, deflected lines concentrate near 
the outer edges of the electrodes ! 43 ' 44 Thus the magni- 
tude of the local external field tangential to the electrode 
surface is substantially larger than Hz- In the analo- 



gous situation of a disc of radius r and thickness a in the 
Meissner state in a uniform transverse field, the field at 
the edge of the disc is a factor of (1 + y t /2r/a) greater^ 

FigurefTOldisplavs the Ic{H) patterns recorded for var- 
ious perpendicular field strengths Hz applied to the same 
initial single- vortex junction state, which was obtained by 
cooling in zero field. An offset in the in-plane field H, due 
to deflection of the applied transverse field Hz through 
the tunnel barrier, has been removed from each Ic(H) 
pattern. The approximately straight lines of the pat- 
terns are associated with the vortex escaping over a po- 
tential barrier due to either the magnetic field H, the 
microshort, or the injectors as indicated in Fig.[7]for zero 
perpendicular field Hz- The measured microshort depin- 
ning current clearly increases with increasing magnitude 
of the circulating supercurrent induced by the perpendic- 
ular field. 

The dipole strength k m depends on the applied trans- 
verse field Hz as k^ = —((Hz + Hzo) where £ is a 
proportionality constant. For small k m or xq, Eq. (|48p 

approximates to U ,J '(xo) = p-sech xq where p, is an ef- 
fective microshort strength that varies linearly with the 
transverse field Hz ■ 



fi = [i + 2((Hz + Hzo) 



(49) 



With this assumption, we fitted the microshort depin- 
ning current to Eq. (|22p and determined the magnitude 
of the background perpendicular field Hzo = 0.011 A/m 
and the constant £ — 8 m/A using the method of least 
squares. An additional fit parameter was the degree of 
misalignment of the solenoid field H. 

The measured microshort depinning current depends 
more strongly on the in-plane field H than expected 
from Eq. (J22) , as is seen from Fig. [7] We believe that this 
was caused by a misalignment of the solenoid field H. As 
is shown in Fig. [101 a misalignment of cither Q\ — 0.057° 
out of (dotted lines) or 62 = 35° within (dashed lines) 
the junction plane fits the observed H dependence of 
the microshort depinning current. In both cases, one 
fit curve is plotted for each of the four values of per- 
pendicular field Hz under which Ic(H) patterns were 
recorded. That the gradient of the measured microshort 
Ic{H) curve depends so sensitively on the degree of mis- 
alignment out of the junction plane results from the 
solenoid field H having a z component when directed 
out of the plane. The effective microshort strength ft is 
then a function of this component H sin 6\ as well as the 
intentionally applied perpendicular field Hz- The larger 
than predicted gradient of the microshort Ic(H) curve 
could also be due to a simultaneous misalignment of the 
solenoid field H both within and out of the junction 
plane. Hence the discrepancy seen in Fig. [7] between the 
measured microshort depinning current and the predic- 
tion of the one-dimensional model, given by Eq. ([^)) . can 
be attributed to a misalignment of the solenoid field H 
in addition to the presence of a small background per- 
pendicular field Hzo = 0.011 A/m. One means of avoid- 
ing interaction between the lithographic microshort and 
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transverse field is described in the next section. 



C. Decoherence sources and improved design 

The vortex qubit has been predicted to have a long 
coherence time for a superconducting qubit, on the or- 
der of tens of microseconds.— This estimate considered 
the effect of quasiparticle dissipation and weakly fluctu- 
ating critical and bias currents in a long linear junction. 
Low frequency critical current noise was found to lowest 
order not to alter the shape of the vortex potential and 
hence not to contribute to decoherence. Since the equi- 
librium density of quasiparticles is exponentially small 
at low temperatures, the coherence time was predicted 
to be limited by low-frequency bias noise. 

Additional factors which could reduce the coherence 
time in our proposed qubit arc increased sensitivity to 
flux noise, the presence of excess quasiparticles, and in- 
teraction with two-level systems in the dielectric. Flux 
fluctuations in the y or z directions modulate the bar- 
rier height of the microshort qubit, causing exponential 
change in the tunneling frequency. Therefore it is im- 
portant to minimize flux noise. Excess quasiparticles 
are generated when the junction switches to the voltage 
state. Their contribution to shot noise in the bias cur- 
rent can be decreased with quasiparticle traps4£ Two- 
level systems are thought to be a major source of de- 
coherence in superconducting qubits4i They originate 
from charge fluctuations and couple to the qubit via the 
electric field ip/tj^- where tj is the tunnel barrier thick- 
ness. In a microshort qubit, two-level systems near the 
microshort may couple to the qubit via the phase oscil- 
lation which results when the vortex oscillates between 
the \L) and \R) states. The amplitude of the phase os- 
cillation is largest at the microshort and decays over a 
few Josephson lengths. The number of two-level systems 
which interact with the microshort qubit could be re- 
duced by employing a tunnel barrier of a-Si:H^ SiN^^i 
or epitaxial A1 2 3 ^ 

Decoherence due to flux noise transverse to the junc- 
tion plane is avoided in a microshort qubit without an 
electrode loop. For example, an improved qubit design 
is a long junction which consists of two linear segments 
separated by an annular segment that contains the litho- 
graphic microshort. At the place where the junction 
ccnterlinc transitions from annular to straight, the field- 
induced potential barrier levels off as seen from Eq. ([6j) 
because, in the straight segments, the component hn 
of the in-plane field transverse to the junction center- 
line is constant. The location of the transition is cho- 
sen so that, at the readout field, the microshort depin- 
ning current 702 is greater than the field depinning cur- 



rent 7c*i. This enables state detection with an RS flip 
flop at one end of the junction-^ provided that the read- 
out bias permits the vortex to overcome the field-induced 
potential barrier but not the microshort-induced barrier. 
At the other end of the junction is a single flux quan- 
tum (SFQ) generator for qubit initialization*^ As well 
as insensitivity to transverse field noise, this design of- 
fers the advantage of faster readout, during which fewer 
excess quasiparticles are generated. Also, the length of 
the junction no longer constrains the radius of the annu- 
lar part; By choosing a shorter radius r, control pulses of 
a given amplitude Ah are achieved with smaller ampli- 
tude pulses of the in-plane field H . Using in-plane field 
pulses large enough that the potential transforms from a 
double to a single well could reduce the sensitivity of the 
coherent oscillation frequency to flux fluctuations in the 
y direction.— 



V. CONCLUSION 

We have studied a vortex qubit based on an annular 
Josephson junction containing a lithographic microshort. 
From the one-dimensional vortex potential, we have de- 
rived the magnetic field dependence of the vortex depin- 
ning current over a microshort-induced potential barrier. 
We have also obtained the frequency of vortex oscilla- 
tion within the well. From this we find that, for a low 
microshort-induced potential barrier, the exponent in the 
coupling Ao between degenerate minima varies with field 
as(l-/i) 3 / 2 . 

The proposed vortex qubit design has been tested 
experimentally in the classical regime. We observed 
bistable vortex states located on either side of the mi- 
croshort. Preparation of the vortex in a given potential 
well was achieved by means of a shaker sequence of bias 
current pulses. We noticed that the depinning current 
from a lithographic microshort can be enhanced by flux 
trapped around both superconducting electrodes as well 
as field applied perpendicular to the junction plane. 
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